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The Levy diffusion processes are a form of non ordinary 
statistical mechanics resting, however, on the conventional 
Markov property. As a consequence of this, their dynamic 
derivation is possible provided that (i) a source of randomness 
is present in the corresponding microscopic dynamics and (ii) 
that the consequent process of memory erasure is properly 
taken into account by the theoretical treatment. 

05.40.+j,05.60.+w 



I. INTRODUCTION 

The theoretical treatment of anomalous diffusion, 
namely diffusion processes either faster or slower than or- 
dinary Brownian diffusion, is an active field of research. 
A well known case of superdiffusion is given by the dif- 
fusion processes of Levy type [jlj . The reader can find an 
exhaustive discussion of the recent literature on this sub- 
ject in excellent review papers However, we would 
like to draw the reader attention also on some key papers 
of the literature on this subject |q-|lS||. We plan to adopt 
a dynamic |l9| , p0| rather than a probabilistic approach 
[|6| — - To make the significance of this purpose more 
transparent it is convenient to compare what we mean 
by dynamic to the conventional probabilistic treatment, 
either resting on the Levy flight or the Levy walk method. 



A. Levy flight and Levy walk 

Both the Levy flight and the Levy walk method are 
based on a totally probabilistic treatment. The Levy 
flight method is based on the assumption that at regu- 
lar time intervals a space transition of arbitrarily large 
intensity might take place. With the Levy walk, on the 
contrary, the jumps over larger distances take place in 
larger times. This property make these processes non- 
Markovian and consequently the derivation of Levy dif- 
fusion more delicate than from within the Levy flight 
perspective. This is easily realized, for instance, by using 
the continuous time random walk formalism PJl^| and 
expressing the time evolution for the probability that the 
particle is at a given space location at a given time by 



means of the equivalent generalized master equation, see 
also [^5[|6). ft is then easily seen that the case where 
the waiting time distribution is characterized by a finite 
time scale yields immediately a Markov process in the 
long-time limit, and the anomalous diffusion properties 
only depend on the long-range nature of the displacement 
per step distribution. In Section II of this paper we shall 
study a physical condition of the same kind. Some special 
attention has to be devoted therefore to the Levy walk 
condition, since it shares with the dynamic approach (see 
Subsection 1C for a more precise definition of this ap- 
proach) a long-time memory, which has to be properly 
erased to establish in the long-time limit the conditions 
for Levy statistics. 



B. From the probabilistic to the dynamic approach 

The dynamic approach, whose precise meaning will be 
discussed in Section 1C, is still somewhat obscure due to 
some conflicting aspects of the recent theoretical deriva- 
tions used to realize this goal. This is so, in spite of the 
fact that the explicit adoption of techniques derived by 
the random walk literature yields a satisfactory deriva- 
tion of the processes of Levy diffusion from within the 
theoretical framework of probabilistic treatments , || , 
0, ||. The purpose of this paper is that of affording a 
unified perspective with no internal contradictions. Here 
we limit ourselves to point out some aspects of the Levy 
walk method which must be retained by the dynamic 
approach. 

First of all, as made clear by the work of Klafter and 
Zumofen 17 and Zumofen ad Klafter jl8| (see also the 
report of Klafter, Zumofen and Schlesinger H), we have 
to point out that the process of Levy diffusion can be de- 
rived from within a dynamic perspective if the so called 
Levy walk view is adopted. This means a trajectory 
moving with constant velocity along a straight line for 
an extended time and from time to time making abrupt 
direction changes. The time of sojourn in one of these 
straight paths is characterized by the probability density 
function 
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where TJ (// — 2) denotes the mean waiting time. The 
renormalization group method, as illustrated by Za~ 
slavsky plj , affords a reliable way of hxing the time 
asymptotic form of (^), and, notably, the power index 
H in terms of the rescaling properties of the fractal re- 
gion at the border between chaotic sea and stability is- 
lands. The theorem of Kac p2j ensures that the first 
moment of tp(t) is finite. This important theorem refers 
to the distribution of the Poincare recurrence times un- 
der the crucial condition that the system under study is 
ergodic. Zaslavsky |Hj noticed that when a stability is- 
land is imbedded within the chaotic sea, the distribution 
of Poincare recurrence times becomes equivalent to the 
distribution of the times of sojourn at the border between 
chaotic sea and stability island. This is so because a tra- 
jectory moving from a given small portion of the chaotic 
sea through a fast diffusion process arrives at the frac- 
tal region and it sticks to it for an extended time before 
returning to the departure region. Thus, in a full accor- 
dance with the Kac theorem, the renormalization group 
theory yields fi > 2, thereby insuring that the first mo- 
ment of this distributions is finite, and consequently that 
also the distribution of the Poincare recurrence times is 
finite. 

We note that all this supports the asymptotic form of 
(0) leaving, though, the impression that its short-time 
structure is arbitrary. It is not so. As proved in a re- 
cent work, gjj , the whole structure of (p is dictated by 
the principle of entropy maximization provided that the 
entropy used is that of Tsallis rather than the con- 
ventional Gibbs-Shannon entropy. Note that the explicit 
form of (|l]) that we are using, is fixed, of course, by both 
the normalization condition and the condition that the 
first moment is finite and that its value is T/ (// — 2). 
In conclusion, the form of Eq.(Q) is a unique analyti- 
cal expression determined by the joint use of dynamics, 
renormalization group technique, and entropy. 

We hope that with no sacrifice of the most important 
ingredients behind the dynamic derivation of Levy dif- 
fusion process, we can restrict our investigation to the 
one-dimensional case. In this condition the role of deter- 
ministic and dynamical generator of the Levy diffusion 
can be properly played by the intermittent map |25[| used 
by Zumofen and Klafter |l^] and by Klafter and Zumofcn 
Jl7| . This is a map with the same algorithmic complexity 
as the Manneville map |^6| , the complexity of which has 
been studied by Gaspard and Wang |^7j by means of the 
Kolmogorov-Sinai entropy. 

In the Hamiltonian model of Zavslasky the derivation 
of the diffusion processes of Levy kind rests on the mi- 
crocanonical conditions. This means that the kinetic 
energy of the flight process is fixed. Consequently, the 
one-dimensional version of this Hamiltonian perspective 
yields the important property that only two velocity 
states exist, one with velocity W and one with velocity 
—W. As a consequence of the one-dimensional assump- 
tion, therefore, we are allowed to use the key relation 
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which is equivalent to 
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Equations (g) and (g) relate to one another the phys- 
ical properties ip(t) and "^(i). The former property , 
ip(t), is the probability density function of sojourn times, 
which, as earlier stressed, has an inverse power law form 
(see Eq.(|l|)); the latter, is the stationary correla- 

tion function of the dichotomous variable £ , playing the 
role of a velocity with only two possible values, W and 
—W. The function ^(i) is determined by the statisti- 
cal properties of the velocity of the paths moving with 
constant velocity and without changing direction. We 
note that Eq. (||) establishes that %p{£) is proportional 
to the second-order time derivative of the function (t) , 
thereby implying, as a consequence of Eq. ([!]), that for 
t — » oo the decay of ^(i) is proportional to 1/t 13 with 
P = fj,— 2. The region of interest for us is that where the 
first moment of ip(t) is finite (so as to fit the Kac theorem) 
and the second moment is divergent so as to prevent the 
system from falling in the attraction basin of the central 
limit theorem. Consequently, we restrict our analysis to 
the interval < < 1. Note that (||) is exact (see p5|) 
if the assumption is made that the time interval between 
the transition from one to the other velocity state is in- 
stantaneous. To help the reader to understand the main 
conclusion of this paper we have also to make another 
preliminary remark. The correlation function is 
a stationary property Jl9[|, implying the existence of an 
invariant, or equilibrium, distribution. A genuinely dy- 
namic approach to the Levy processes, consistent with 
the ergodic assumption behind the Kac theorem, implies 
that this equilibrium distribution is established by a sin- 
gle trajectory, provided that this trajectory runs for an 
unlimited amount of time. The lack of a finite micro- 
scopic time scale makes this condition difficult to realize 
in practice, and it is probably one of the sources of the 
conflicting views that will be discussed in this paper. 



C. The dynamic approach 

The general program of the dynamic approach to sta- 
tistical mechanics is illustrated in a series of recent papers 
p8| . We are very close to the program of Ref. . 
The ambitious purpose of these authors is to derive an 
important equation as the Fokker-Planck equation with- 
out using any statistical assumption whatsoever, so as 
to reverse the ordinary path from thermodynamics to 
statistical mechanics. In other words, the path to fol- 
lows moves from dynamics and reaches the level of sta- 
tistical mechanics using only deterministic randomness 
with no recourse to thermodynamics, this being the last 
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step, stemming from the dynamically generated statisti- 
cal equilibrium distributions. 

The authors of p3| adopted the same perspective to 
move from dynamics to Levy statistics. The authors of 
Rcf. [M found that the density distribution a(x,t) of the 
variable x driven by a process described by Eq.(g) obeys 
the equation of motion : 

Jl t) = U 2 > J* $ c (f) a(x, t - f ) dt'. (4) 

Within the context of a dynamic approach to the Levy 
processes, this equation should be given a special atten- 
tion, since no explicit use of probabilistic arguments was 
made to derive it [^0|. However, no general solution of 
it is available, and the emergence of the Levy diffusion 
out of it rests on an approximation which has been ques- 
tioned Jp3| , (36[ : Different approximations to the solution 
of Eq.(plead to different statistical processes. The in- 
terested reader is referred to the work of Ref. [^0| for 
the the derivation of Eq.(|J). Here we limit ourselves to 
noticing that this equation is exact under the condition 
that the velocity variable is dichotomous and the initial 
distribution is a Dirac delta centered at x = 0. Thus, 
there is an intimate relation between Eq.(|J) and Eq.(|j). 



D. Purpose and outline of the paper 

The main purpose of this paper is that of stressing 
that the intimate connection between Eq. (|2|) and Eq. (Q) 
does not leave room for a solution implying a memory 
infinitely extended in time. In other words, we want to 
prove that the adoption of a Markovian perspective, al- 
though apparently incompatible with the time convolu- 
tion of Eq. (Q), is dictated by the steady action of the 
randomness corresponding to the transition from the one 
to the other state of the variable velocity. Consequently 
the Markov structure emerging from Eq. (^J) according to 
the prescriptions of p(| cannot be misled as an undue 
approximation. This is rather an ingenuous way of estab- 
lishing a physical condition fitting the result of an earlier, 
and crucial, research work. Gaspard and Wang Q prove 
that in the long-time limit the Manneville map becomes 
equivalent to the Bcrnouilli shift map. This is a conse- 
quence of the fact that the repeated exit and re-entering 
into the laminar region results in a memory erasure. As 
we shall see, this is the main reason why the final state 
is as Markov as a genuine Levy process must be. 

In other words, we plan to make randomness emerge 
from the dynamic approach, resting on Eq.(^), so as to 
render this dynamic approach equivalent to the Levy 
walk perspective. The purpose of this paper is to show 
that the Markov property necessary to derive the process 
of Levy diffusion is not arbitrary, but rather corresponds 
intimately to the nature of the dynamic process resulting 
in Eq.(|4|). This is so because Eq.([I|) implies the dichoto- 
mous nature of the variable £. The condition where for 



t — > oo the decay of 3>f(i) is proportional to l/t& with 
[3 = /i — 2 means that the process is deterministic for 
an overwhelming amount of time. There exists an inti- 
mate equivalence between Eq.(^|), in this physical condi- 
tion, and intermittent maps. Randomness shows up only 
when the trajectory reaches the border between laminar 
and chaotic region j26). At this crucial stage there are 
only two possible events, either a jump from the original 
into the other laminar region, corresponding to a distinct 
velocity state, or the jump back to the original laminar 
region, namely, the original velocity state [ ji8| . At this 
stage dynamics are essentially indistinguishable from the 
time evolution of the Bcrnouilli shift map, whose con- 
nection with thermodynamics and statistical mechanics 
has been recently clearly illustrated by Zaslavsky 
This means that randomness is a rare event and it is 
in fact the reason why we have adopted the concept of 
sporadic randomness. The main purpose of this paper 
is that of making a choice between two distinct ways of 
solving Eq.(|j), based on the criterion that the right so- 
lution must reflect this sporadic randomness. 

The outline of this paper is as follows. In Section II we 
review the arguments used in an earlier paper, |35| ], to de- 
rive a process of Levy diffusion by means of a generalized 
master equation. In Section III, using the calculation il- 
lustrated in the Appendix, we show that the same result 
is derived from a master equation which looks like the 
Markov approximation of that of Section II. In Section 
IV we review, in the light of the perspective established 
in this paper, the method used in Ref. to assign to 
Eq(.^j) a Markov structure. Some final conclusions are 
made in Section V. 



II. THE GENERALIZED MASTER EQUATION 

The first step of our approach rests on the use of the 
generalized master equation of Rcf. |57|. This equation 
reads: 



d_ 

at 



o~(x, t) 



dt' 



n(x-x',t- t')o{x', t') dx' , (5) 



k(x, t) = tt(x, t) — 8{x) I dx'ir(x',t) (6) 



where 



and ir(x,t) denotes the probability for the particle to 
make a jump by a distance x at time t. This equation is 
very general and is expected to be compatible with the 
description of highly non-Markovian processes like that 
corresponding to Eq. (Q) with the non integrable correla- 
tion function generated by Eqs.(|f|) and (||). The intimate 
connection between these two equations will be discussed 
in Section IV. Here we limit ourselves to stressing that 
the asymptotic regime of a(x,t), as given by Eq.rtq), can 
be studied without making explicitly the Markov approx- 
imation. In fact, using the property that this equation 



3 



is convoluted in both space and time variables, we get 
for the Fourier-Laplace transform of a(x,t), denoted by 
a(k,s), the following expression: 



H 



a(k, s) = 



1 



s — k(k, s) 



(7) 



where, of course, £(fc, s) denotes the Fourier-Laplace 
transform of /v(ir,i). As pointed out in Section I, the 
dynamic approach to Levy statistics that we are consid- 
ering is related to the Levy walk condition. This means 
that a transition of length \x\ implies a time t = \x\/W. 
In conclusion we are forced to make the following choice 
for tt(x, t), with t > 0: 



n(x,t) = tp(t)S{\x\ - Wt). 



(8) 



The authors of Ref. |33| studied the asymptotic regime of 
Eq. (Q) , supplemented by Eq. (||) , searching for the rescal- 
ing condition 

socfc Q , (9) 

with a > 1. This is a reasonable assumption, since in 
the asymptotic limit the second moment is known [ p3| to 
yield 



2H 



with: 



H = 1-/3/2 



(10) 



(11) 



On the other hand the rescaling of Eq. (|]) suggests that 
the condition a — 1/H might apply, thereby resulting in 
the property a > 1 , which is essential for the calculations 
aiming at establishing the exact dependence of a on /?. 

With straightforward calculations it is shown |33| that 
in the asymptotic limit Eq.(^|) and Eq.(||) yield: 



a = p+l 



and 



with 



a(k, s) 



b\k\ a ' 



(WT) a an 
b = T(l-a)^-(a-l)coa(^). 



(12) 



(13) 



(14) 



Note that this is the Laplace transform of the following 
equation of motion: 



J^o-(/M) = -b\k\ a a(k,t). 



(15) 



This means that the asymptotic regime of the generalized 
master equation of Eq. ( 5 ) is a process of diffusion with a 
genuinely Levy nature [2|. 

We note that Eq.(Q2h means the rescaling 



1+0 



(16) 



which is different from that of Eq. (|l 1[) . The difference 
between the two rescalings is a fact of crucial importance 
deserving proper comments. We note that the rescaling 
of Eq. (^TJ) is somewhat ambiguous since it refers to the 
dynamics of Eq. (0) . As pointed out in Ref. [|TJ , as well 
as in the earlier work of Zumofen and Klafter jlj| and 
Klafter and Zumofen j[7) , the diffusion process described 
by Eq.(^) consists of a central part and a propagation 
front signalled by two sharp peaks. At time t a particle 
leaving the origin x = at t — cannot be found at a dis- 
tance from the origin larger than Wt. This has the effect 
of producing an accumulation of particles at the front it- 
self of the diffusion process, namely at x — +Wt. This 
is the origin of the two ballistic peaks of the propagation 
front. At earlier times the initial distribution, concen- 
trated at x = 0, splits into these two ballistic peaks and 
the region between the two peaks is empty. Due to the 
effect of sporadic randomness, some trajectories leave the 
propagation front and the population of the central part 
steadily increases in time, while the peak intensity, pro- 
portional to the correlation function <£>£ slowly decreases. 
Note that this means that the diffusion process cannot be 
described by a single rescaling. The peaks of the propaga- 
tion front rescale with H = 1, a fact implying a diffusion 
faster than that predicted by the rescaling of Eq.(^Tj). 
The rescaling of the central part is properl y e xpressed 
by Eq. (|l6|) . The calculations leading to Eq.(^5|) refer to 
a physical condition where the intensity of the ballistic 
peaks is negligible, so that the rescaling of Eq.(|l6|) only 
reflects the diffusion properties of the distribution central 
part. On the contrary, the rescaling of Eq.(|Tl|) is a sort 
of balance between the fast rescaling of the propagation 
front and the rescaling of the central part of the distribu- 
tion a{t) , which is in fact slower than rescaling of Eq. ([16]) . 
In conclusion, the discrepancy between the rescaling of 
Eq.(|l"T|) and the rescaling of Eq.(16) reflects the fact the 
derivation of the Levy process from within the Levy walk 
perspective rests on a deep conflict between the dynamic 
properties still present within the Levy walk perspective 
and a merely probabilistic treatment. 

We want to make a further remark, concerning the 
derivation of a Levy process from the generalized master 
equation of Eq. (||) . We note that the Levy processes are 
a form of Markov statistics [Q thereby implying that the 
asymptotic regime of the dynamic process described by 
Eq. (|) involves a process of memory erasure. This pa- 
per is devoted to proving that Eq. (|J) , in spite of its non- 
Markovian structure, yields an asymptotic regime whose 
statistics are determined by this process of memory era- 



4 



III. THE MARKOV MASTER EQUATION 

As pointed out at the end of Section II, the processes 
of Levy diffusion are Markovian. Therefore, it is conve- 
nient to discuss their derivation from a Markov master 
equation. This discussion will serve the useful purpose of 
proving that the master equation is, in a sense, a bridge 
between the dynamic treatment and the diffusion regime. 
The latter is the subject of thermodynamical arguments 
and the former rests, in the theoretical picture adopted 
in this paper, on classical mechanics. Thus, the master 
equation can also be regarded as an important bridge 
between mechanics and thermodynamics. In the contin- 
uous representation the master equation reads, see [Q, 
§1 and M: 



d_ 

dt 



cr(x, t) 



K(x - x')a{x' ,t)dx' 



where: 



(17) 



(18) 



As far as the transition probability H(x) is concerned, 
we adopt the result of the entropic analysis of Ref. [^3) . 
Thus we write: 



K(x) = n(x) - 6(x) / IL(x')dx 



y J TW V \W ) {TW + \ x \y 



(19) 



with 2 < /i < 3. In fact, the theoretical work of Ref. [23] 
proves that the maximization of the Tsallis entropy [24] 
with entropic index q = 1 + 1/fx yields for II an inverse 
power law form with index fj,, provided that a transition 
of length |ir| is supposed to be related to the time t by 
\x\ = Wt.The adoption of this entropic argument changes 
(O) into: 



+ 00 



\x-x'\ 2 +P 



o~{x' , t)dx' 



(22) 



and, in this sense, coincides with the expression found by 
Seshadri and West M. Note that keeping T > makes 
it possible for Eq.(j20|) to cross the critical condition (3 = 
1 without meeting the divergence corresponding to the 
Levy prescription, namely the divergence of b of Eq.(|l4|) 
at a = 2. The result of Eq.(|2(i|) can be used also to study 
the region (3 > 1 corresponding to the attraction basin of 
the ordinary central limit theorem. 

It is important to observe that the Markovian master 
equation here under study can be derived from the gen- 
eralized master equation of Eq.(||) by using the Markov 
condition: 



1%) 



n{x,t')dt'. 



(23) 



This is an axpect of crucial importance. In fact Eq.(|2E 
yields Eq. ( |l9| ) , showing that the Markov property makes 
it possible to establish an important connection between 
dynamics and thermodynamics. The Markov perspec- 
tive adopted in this Section is essential to establish the 
key connection between the structure of Eq.(|l9|), result- 
ing from the adoption of entropic arguments, and the 
structure of Eq.(|l]), generated by the adoption of the 
renormalization group arguments, which, in turn, reflect 
genuinely Hamiltonian properties pl| . Note that we re- 
fer to Eq.(p3|) as a Markov condition rather than as a 
Markov approximation. This is so because the term "ap- 
proximation" suggests a given departure from the exact 
solution or, in other words, an error whose intensity must 
be defined. We see, on the contrary, that the asymptotic 
regime of Eq.(fj"o|) coincides with that of Eq.(g), if, the 
latter equation is supplemented by the crucial condition 
of Eq.(||), mirroring the dynamics illustrated in the In- 
troduction. 



d_ 

di 



a(x, t) = 



( P -i)(Twy 



T 



a(x' , t) dx' 



{TW 



x'\) 



2+/3 



dx' 



{TW - 



\) 2+f} 



a{x, t) 



(20) 



where fi — 2 + with < ft < 1. As shown in Appendix, 
the Fourier transform analysis of ( pfj| ) proves that in the 
asymptotic regime {kWT -C 1) this equation is equiva- 
lent to the process of Levy diffusion of Eq. (jl5|) where the 
parameter b reads: 



, „ / 7TQ! \ , (WT) a 

6 = 2cos(— )r(l- a )^J- 



(21) 



In other words, we derive a result equivalent to that of 
Section II. Note that in the limiting case T = the term 
within the square brackets of Eq. (EG) becomes the regu- 
larized [B9I1 form of: 



IV. THE EXACT DIFFUSION EQUATION AND 
THE MARKOV REGIME 

Is there a connection between Eq.(^) and Eq.(po|)? It is 
evident that this connection would be established by the 
exact solution of (0), if this yielded, in the asymptotic 
time limit, a diffusion process of Levy kind. It has to 
pointed out, however, that finding an exact solution of 
(4) is not easy. In literature we find only solutions of 
(4) based on approximations p0p^ , |3^ ]. In Ref. [M a 
solution of Eq.([|) was found, with a Markov character, 
and corresponding to Eq.(20) with only the first of the 
two terms between the square brackets on the r.h.s of 
this equation. In ref. |35) a non-Markov solution has 
been discussed, which has been later judged to be the 
correct solution fl36| . 

The purpose of this Section is that of discussing these 
two proposed solutions in the light of the sporadic ran- 
domness illustrated in Section I. In the work of |3q] it 
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was argued that an analytical solution of Eq. Q4p can be 
found, based on the fact that a fractional derivative in 
time emerges from the r.h.s of (||) if the correlation func- 
tion $£(t), which should fulfill the normalization condi- 
tion $£(0) = 1, is replaced by a function like const /t". 
This means a function with the same long-time property 
as the original correlation function, breaking however the 
normalization condition at t = 0. This approximation re- 
sults in very appealing mathematical properties. In fact, 
it has the nice effect of resulting in a process with in- 
finite memory, and in an analytical expression for the 
effects that this infinitely extended memory has on diffu- 
sion. The adoption of this approximation yields p6| the 
rescaling of Eq. (|TT|) . It is straightforward to prove that 
this rescaling can be obtained from the Fourier-Laplace 
transform of Eq. (H) 



a(k, s) 



<e> 

+ <£ i (s)k 2 



(24) 



where ${(s) is the Laplace transform of the correlation 
function ${(t). To derive the rescaling of Eq. (|ll| ) 
we have to replace in Eq. (HJ) $g(s) with s 13 " 1 (see 
p8| ). This rescaling, however, conflicts with the numer- 
ical observation of Ref . |^o| , which results in a different 
rescaling, corresponding to that of Eq. (|l6|) . The dis- 
cussion of Section II sheds light on the origin of this 
rescaling, different from that of Eq. (p"l|). It seems 
to be evident to us that the study of the asymptotic 
properties of Eq. ( p4| ) resting on the limiting condition 
lim s _ >o $j (s) = const s' 3-1 , loses any dependence on the 
key parameter T, and with it, on the fact that there ex- 
ists a propagation front moving with finite velocity. This 
explains why the same method of time asymptotic anal- 
ysis applied to Eq. (0) yields the correct rescaling. This 
is so because in this case the Levy walk nature of the 
process under study is retained by the kernel n{x,t) due 
to the wise choice made for tt(x, t) in Eq. (||). 

In conclusion, we are convinced that the solution im- 
plying the existence of an infinitely extended memory 
conflicts with the numerical treatment of the diffusion 
process resulting from the fluctuations of the dichoto- 
mous variable £ with a non integrable correlation function 
$f(t). This is so because the steady action of sporadic 
randomness has the effect of producing a Markov statis- 
tics, although this occurs in the long-time limit. How to 
make the rescaling of the central part of the distribution 
a(t) to become compatible with the effect of sporadic ran- 
domness and with the predictions of the Levy-Gnedenko 
theorem in the long-time limit? The most direct way 
to realize the correct rescaling of the central part of the 
distribution and to make the Markov property emerge is 
that of assuming that: 



1 f°° 

a(x, t-t') = - I 8{W t' - \ x-x'\) a(x', t) dx' 

* J-oo 



As earlier pointed out, if we apply this condition to the 
term on the r.h.s. of Eq.(^), we obtain an equation of 
motion identical to one that would result from Eq. (pp|) 
cancelling the second of the two terms within the square 
brackets of this equation. This interesting result implies 
some algebra based on the method of integration by parts 
and the properties of the delta of Dirac. More details on 
this calculation can be found in Ref. p(i|| . 

It is important to point out the physical meaning of 
the constraint of Eq. (|25|) . This means that we imag- 
ine a condition still unaffected by randomness, since this 
constraint would be rigorously valid only in the case of 
merely ballistic motion. Yet, the effect of replacing ( j25| ) 
into the r.h.s of (|4|) is that of producing a Markov struc- 
ture as an effect of carrying out the integration on t' in 
the time convoluted form of Eq.(^). We are convinced 
that the emergence of this Markovian structure is not an 
artefact of the approximation of Eq. fl25| ) . The error asso- 
ciated to this approximation is not the emergence of the 
Markov structure. This error is totally different in na- 
ture, and can be easily evaluated. In fact, as repeatedly 
pointed out earlier, this approximation has the effect of 
resulting only in the first term on the r.h.s of (|20|). The 
error associated to this approximation is signalled by the 
breaking of the norm conservation. It is evident in fact 
that the condition: 



cr(x, t) dx = 1 



(26) 



is fulfilled by both (|j) and ( pcf ) and that in the latter case 
this is the consequence of the wise structure of the mas- 
ter equation of (p"7|). This is the reason why we look at 
the master equation of ( po| ) as a natural bridge to cross 
when moving from the dynamical perspective of (Q) to 
the final regime of Levy kind. As noticed in Section III, 
the emergence of this final condition from (|2^) is made by 
using the mathematical approach of the Appendix which 
proves that the Fourier transform of ( p0|) yields the gen- 
uinely Levy process of diffusion of Eq.([15|). We note that 
in the earlier work 20 . no proper attention was de- 
voted to the crucial fact that the second term within the 



(25) 



square bracket of Eq. (|20| ) is essential for a proper deriva- 
tion of the Levy processes, or of the equivalent fractional 
derivative. 



V. CONCLUDING REMARKS 

In conclusion, we have provided a convincing demon- 
stration of how to derive Levy processes from within a 
dynamic approach. As pointed out in an earlier work 
p0| the Levy process corresponds to a form of fractional 
calculus which can be regarded as a form of macroscopic 
manifestation of microscopic randomness. However, it 
seems to us that fractional derivatives in time have a dif- 
ferent meaning from fractional derivatives in space. The 
choice made by the authors of Ref. |3q] has the effect 
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of relating the solution of Eq. (Q) to a form of fractional 
derivative in time. It seems to us that this choice cor- 
responds to a case where the decay of the correlation 
function is not originated by the sporadic action of ran- 
domness on a single trajectory. In this latter case, as 
shown in this paper, to obtain the correct result we are 
forced to make the Markov property emerge, and this is 
realized by adopting the trick of Eq. (p5|) . The assump- 
tion made by the authors of Ref. |3^], on the contrary, 
seems to imply that the decay of the correlation function 
is originated by a statistical distribution over a range of 
initial conditions. Whether or not this is compatible with 
the dichotomous nature of the diffusion generator is not 
quite clear to us. 

Should the proof be given that the relaxation of the 
correlation function $j can also be determined by statis- 
tics as well as by dynamics, with no conflict with the 
dichotomous property behind Eq. (Q), we would reach 
the impressive conclusion that this exact equation admits 
two distinct classes of solution, determined by the extra 
information about the physical origin of the relaxation 
process. We are inclined to believe that the relaxation of 
a dichotomous variable can only be compatible with the 
action of a sporadic randomness. Thus the emergency 
of the Markov property at extremely long time is a fair 
reflection of the dynamical processes generating the Levy 
diffusion processes, not to speak of the fact that the corre- 
sponding rescaling of the central part of the distribution 
fits the results of the numerical calculations |^(| . 

It is also important to stress that the repeated action 
of randomness is subtly related to the possibility of es- 
tablishing a connection between mechanics and thermo- 
dynamics also in the case, studied in this paper, of ap- 
parently infinite memory. The random seed is given by 
the fact that the duration of the times of sojourn in the 
laminar region cannot be predicted. This is so because 
of the random injection into the laminar region from the 
chaotic part of the map j25j . This is the reason why we 
have to introduce probabilistic arguments within the dy- 
namical picture of the process under study. This is also 
the reason why, as shown by [^3| , the shape of the density 
distribution of Eq. ([!]) can be predicted by using entropic 
arguments, provided that the non-extensive form of en- 
tropy advocated by Tsallis (24| is used. In other words, 
both the adoption of entropic arguments and the birth 
of Levy statistics rest on the emergence of probabilistic 
aspects generated by a sporadic form of randomness. 

Finally, we are left with the intriguing issue of the de- 
pendence of all these properties on the space dimensions. 
The treatment of this paper has been confined to the one- 
dimensional case, where the crucial action of the stabil- 
ity islands pointed out by the theoretical analysis of Za- 
slavski [^lj is correctly mirrored by the dichotomous na- 
ture of the fluctuating variable £. We are convinced that 
the essence of the present treatment can be extended to 
the multidimcnsiuonal case, including the more realistic 
three dimensional case. However, we have to recognize 
that this extension is not a trivial matter and that further 



research work has to be done to settle the technical prob- 
lems triggered by the two- and three-dimensional case. 

APPENDIX 

We define the function 



Tin + 1 - a) 

n=0 v ; 



(27) 



This function is a generalization of the exponential func- 
tion. In turn, this generalized exponential is used to de- 
rive the following form of generalized trigonometric func- 
tion: 



and 



cos a x 



2i 



+ E~ 



The corresponding expansion in a power series are 
x™- Q sin[(n - a)ix/2] 



sin 



n=0 



r(n + 1 - a) 



and 



cos a x 



E 



x"- a cos[(n - q)tt/2] 
r(n + 1 - a) ' 



(28) 



(29) 



(30) 



(31) 



Note that the functions defined in Eq.(g7|) fulfill the im- 
portant relation 



J x a e x dx = T(l + a)e- x E x _ x _ a . 



(32) 



We are now ready to address the problem raised by 
Eq. ( pp[ ) , namely,the evaluation of the Fourier transform 
of the function 1/(TW + \ x — x' |) 2+/3 . For this purpose 
we can use Eq. (|32|) and we find: 



+ OC 



(a+\x\) 



On the other hand 
f(k) = 2 



2+/3 



dx = 2Re 



/(*)• 



zkx dx 



{a+\x\y+ a 



(33) 



cos fca 



(a + x) 1 



-dx 



2T(-a) \k\ a [sin[-(l + a) + \ka\] 



-sin a [—(1 + a) + | ka | 



(34) 



where a = WT and 
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sin Q [— (1 + a) + | ka | ] = cos [— (a + 1) ] sin a \ ka | 

+sin[— (a + l)]cos a \ ka |. (35) 

We are interested in the case ka — > 0. Thus we use 
(|30|), ( |3l| ) and (|35|) to evaluate the Fourier transform of 
interest keeping only the leading vanishing and diverging 
terms. We thus obtain: 



+ °° exp(ikx) 
ax 



, , * 2|/fc| Q r(-a)[cos(^) 



•sin(^-)|fca| 



r(l-a) r(3-a) 
In conclusion, the Fourier transform of (^0|) yields 
d 



dt 



a(k,t) = c[f(k) - f(Q)]a(k,t), 



where f(k) is given by Eq.(BJ) and 



(36) 



(37) 



(38) 



where /i — 1 = a. After some algebra we find that /(0) 
2/(aa a ) end the final coefficient in front of a(k,t) is: 



,ira. , . a I ka \ c 
2cos(— )r(-a)^ r l 



, 7ra . „ „ . a , 
-2cos(— )r(l-a)— | fc 



(39) 



As a final result we get Eq. (|l5|) . 
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